The Critical Exponent of the Fractional Langevin Equation is ac ~ 0.402 
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We investigate the dynamical phase diagram of the fractional Langevin equation and show that 
critical exponents mark dynamical transitions in the behavior of the system. For a free and harmon- 
ically bound particle the critical exponent = 0.402 ±0.002 marks a transition to a non-monotonic 
under-damped phase. The critical exponent q_r = 0.441... marks a transition to a resonance phase, 
when an external oscillating field drives the system. Physically, we explain these behaviors using 
a cage effect, where the medium induces an elastic type of friction. Phase diagrams describing the 
under-damped, the over-damped and critical frequencies of the fractional oscillator, recently used 
to model single protein experiments, show behaviors vastly different from normal. 
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Anomalous subdifFusion (x^) ^ with < a < 1 is 
found in diverse physical systems ranging from charge 
transport in disordered semiconductors, and quantum 
dots, to relaxation dynamics of proteins and diffusion of 
mRNA in live E-coli cells to name a few examples [l|, 0, S, 
S, S S S S] • On the stochastic level, anomalous diffusion 
and relaxation is modeled using the fractional Fokker- 

lOj, a,nd the frac- 
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Planck-Kramers equations [3, |S , | 
tional Langevin equation [^.T^.llll. 1 
For example recent single molecule experiments on pro- 
tein dynamics 0, [IH revealed anomalous Mittag-Leffler 
relaxation. The recorded anomalous behavior: the dis- 
tance between an electron donor and acceptor pair within 
the protein, bounded by a harmonic force, was succcss- 
fuUy modeled by the fractional Langevin equation (FLE) 

For the standard Langevin equation with white noise, 
the well known under-damped and over-damped phases 
and the critical frequency yield a dynamical phase dia- 
gram of the motion. In this manuscript fractional over 
and under-damped behaviors are investigated, which ex- 
hibit rich behaviors vastly different from the normal case. 
The most surprising result is the appearance of a criti- 
cal exponent for the FLE - etc ~ 0.402. For a < ac the 
over-damped behavior totally disappears, and the decay 
to equilibrium is never monotonic. We will interpret this 
behavior in terms of a cage effect, and show that simi- 
lar critical behaviors are found also for a particle free of 
deterministic external forces, and for the response of the 
system to an external time dependent field. Thus criti- 
cal exponents control fractional dynamics. Contrary to 
expectation, our work shows that the transition between 
normal diffusion and relaxation a = 1 to the strongly 
anomalous case a — > 0, is not a smooth transition. 

The FLE is a generalized Langevin e quat ion [2, 
with a power-law memory kernel 0, 0, flllll^- 0, 0, 
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external force field, < a < 1 is the fractional exponent 
and ^(t) is a stationary, fractional Gaussian noise 0, \^ 
satisfying the fluctuation-dissipation relation Q 



(2) 



The FLE Eq. (O can be derived from the Kac- 
Zwanzig model of a Brownian particle coupled to an Har- 
monic bath [17|. And recently the non-trivial fractional 
Kramers escape problem for the FLE was solved [ist . 
Another convenient way to write Eq. ([T]) is 
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^F{x)~-fT{l~a) — +m 
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where the fractional derivative [H is defined in Caputo 



sense 



'fit) _ 



^jl[t-t)-'^{df{t)/dt)dt. 



In re- 
cent years fractional calculus was shown to describe many 
Physical systems |8| . 

The normalized correlation function 
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for a harmonic force field F{x) = —mw'^x, was used to 
describe single protein dynamics [1, [lH , so we will first 
focus on this experimentally measured quantity. Ther- 
mal initial conditions are used {^{t)x{0)) = 0, (a;(0)^) = 
kbT/muo^ and (a;(0)w(0)) = 0. We assume that a — p/q, 
where q > p > Q are integers and p/q is irreducible (i.e 
not equal to some other l/n where I < p and n < q are 
integers). Using the convolution theorem and Eq. 
the Laplace transform of Cx{t) is 
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where 7 = ^7r(l — a). We rewrite Eq. ^ as 
(s + 7S?"^)Q(s) 



where 7 > is a generalized friction constant, F{x) is an 



Pis) 



(6) 
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with 



and 



p{s) = {s^ + oj^Y + {-ly- 



(7) 
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By finding the 2q zeros of P(afe) = 0, a^. with fc = f , ...2q, 
and using analysis of poles in the complex plane we 
invert Eq. ^ and obtain an explicit analytical so- 
lution for Cx{t) [131 • Defining the constants Ak — 
1/ dP{s) / ds\s=ak and Brnj by the expansion Q{s){s 



') = E'n^o El=d ^mjs— ^"/«, the solution is 
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liakt), (9) 



where is a generalized Mittag-Lefflcr function l]. 

The solution Eq. ((9]) is valid provided that all the a/jS 
are distinct, otherwise a critical behavior is found which 
is soon discussed. We see that finding the solution to the 
problem is equivalent to finding the zeros of the polyno- 
mial P{s), once these zeros are found one can investigate 
and plot the solution with a program like Mathematica. 
For example for a = 1/2 Eq. ^ reads 



Cx{t) = ELi ((-7' + ^^cik + 4) Ate- 
+-1uj'^Akt^E^^^{akt)) . 



(10) 



where at are non-identical solutions of P{s) — (s^ -I- 
w^f - 7^5 ^ 0. 

The t — > 00 asymptotic behavior is found using the 
exact solution Eq. ^ 



C^{t) 



a] 
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or with Tauberian theorems [2l|. Hence for very long 
times Cx{t) always decays monotonically. However for 
shorter times the picture is very different. In Fig. [U 
the analytical solution for Cx{t) with a = 1/2 is plotted 
for various uj using natural units 7=1. Three types 
of behaviors exist (i) Monotonic decay of the solution - 
Fig. [T]^a) . (ii) Non-monotonic decay of the solution with 
no zero crossing Cx{t) > Fig. [Tfb). (iii) When the 
frequency uj becomes large Cx (t) exhibits non-monotonic 
decay with zero-crossings Fig. [TJc). Similar types of 
behaviors are found also for other parameter set. 

Our goal now is to clarify the phase diagram of the 
dynamics of the fractional oscillator. For a — 1, the reg- 
ular damped oscillator, Uc = 7/2 defines the critical fre- 
quency which marks the transition between over-damped 
and under-damped motion. For the fractional oscillator 
< a < 1 the identification of a single critical frequency 
is not possible, and we use three definitions for critical 




FIG. 1: The correlation function Cx{t) versus t with a = 
1/2 and natural units 7 = 1. Three types of solutions are 
presented (a) uj — 0.3 and the correlation function decays 
monotonically. (b)The critical frequency ui — ujz ~ 1.053 
where the transition between motion with and without zero 
crossing, namely Cx{t) = for a single point in time, (c) For 
larger frequency a; = 3 we observe non-monotonic decay with 
zero crossing. 



frequencies where the dynamics of the system undergoes 
a transition in its behavior. For the normal case a = 1, 
the critical frequency lOc is found when two zeros of P{s) 
coincide, i.e. an appearance of a pole of a second order 
for Cx{s). Similarly for the < a < 1 case we define 
LJc as a critical point for which not all the zeros of P{s) 
are distinct and as a result for such lOc Eq. Q is not 
valid. The second approach is to take the minimal fre- 
quency LUz at which the solution Cx{t) crosses the zero 
line. Finally we will investigate the minimal frequency 
LJm at which Cx (t) is no longer a monotonically decaying 
function. Only for the regular damped oscillator a — I, 
we have ujc = = ^m- 

Eqs. (|5l9p are used now to investigate the phase di- 
agram of the motion. By its definition, at the critical 
frequency Wc a pole of second order in Eq. ([5]) is found. 
It is easy to see that for such loc two conditions must 
be satisfied, P{s) ~ and dP{s)/ds — 0. These two 
conditions lead to the following relation for lOc 
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(2 — a)a2-Q-)/i 



(12) 



which is valid only for even q or even [q + p) where a — 
p/q. For odd q and even p, the critical point lUc doesn't 
exist at all. Similar to Eq. (|12p we find from dimensional 
analysis 



(13) 



where Kzict) and Kmia) depend only on a. By investi- 
gating the exact analytical solution Eq. ^ for various 
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FIG. 2: The phase diagram of the fractional oscillator. Phase 

(a) monotonic decay of the correlation function Cx{t), phase 

(b) non-monotonic decay without zero-crossing and (c) oscil- 
lations with zero crossings. The boundary between (b) and 

(c) is ujz = Kz(q) (solid line -I- squares), the boundary be- 
tween (a) and (b) is ujm ~ f^m{a) (solid line -|- circles). For 
a < ac — 0.402, the phase of monotonic decay disappears, 
namely we do not find over damped behavior. The dotted 
curve is the critical line u>c given by Eq. (|12|l . All the curves 
are calculated for 7 = 1. For a = 1, ujc = uJz = ijJm = 7/2. 



a and 7 = 1 we obtain the functions Kz(a) and Km{oi)- 
The resulting phase diagram of the fractional oscillator 
is presented in Fig. O 

From the phase diagram Fig. [2] a very interesting re- 
sult emerges. For a < ac ~ 0.4 the solution is never 
decaying monotonically. For such a an oscillatory behav- 
ior is always found even if the frequency of the binding 
Harmonic field lu 0. In this sense the solution is al- 
ways under-damped. Thus we find a critical exponent ac 
which marks a pronounced transition in the behaviors of 
the solutions of the fractional oscillator. We also note 
that for a near ac we find Km{a) oc (a — ac)^ , which 
describes the boundary between the monotonic and non- 
monotonic phases. Note that the phase diagram Fig. [2] 
also exhibits some expected behaviors: as we increase 
uj we find a critical line above which the solutions are 
non-monotonic and exhibit zero crossing (similar to usual 
under-damped behavior) and for a = 1 the three critical 
frequencies lOc, oJz and LOm are identical. 

A physical explanation for the phase diagram is based 
on the cage effect. For small a the friction force induced 
by the medium is not just slowing down the particle but 
also causing the particle a rattling motion. To see this 
consider the FLE Eq. ^ in the Hmit a — > 



mx 



(14) 



where xq is the initial condition. Eq. (fHl) describes 
harmonic motion and the "friction" 7 in this a — s- 
limit yields an elastic harmonic force. In this sense the 
medium is binding the particle preventing diffusion but 



forcing oscillations. In the opposite limit of a — s- 1 

mx + m'yx + muj'^x ^(t) when a ^1, (15) 

the usual damped oscillator with noise is found. So from 
Eq. for a ^ an oscillating behavior is expected, 
even when — s- which can be explained by the rat- 
tling motion of a particle in the cage formed by the sur- 
rounding particles. This behavior manifests itself in the 
non-monotonic oscillating solution we have found in our 
phase diagram Fig. [5] when a — *■ 0. Our finding that ac 
marks a non smooth transition between normal friction 
a — » 1 and elastic friction a ^ is certainly a surpris- 
ing result which could not be anticipated without our 
mathematical analysis. 

To find accurate values of Oc we used also a method 
based on Bernstein theorem [2l|. According to the the- 
orem if and only if a function f(t) is positive then for 
any integer n, < (— l)"(d"/(s)/ds"'), where /(s) is the 
Laplace pair of f{t) . To see when Cx (t) is decaying mono- 
tonically, we investigate the positivity of its derivative. 
Using Bernstein's theorem inspecting the n — 150 or- 
der derivatives gives the exact lower bound ac > 0.394, 
while the method using exact solution for Cx{t) gives 
ac = 0.402 ± 0.002. Further we can show that the criti- 
cal exponent is not equal to 0.4, and hence it seems a non 
trivial number. We now show that the critical exponent 
ac is important for other physical quantities. 

The mean square displacement {x'^{t)) for the force 
free particle F{x) — is now investigated. We con- 
sider averages over the noise and thermal initial con- 
ditions (w^)|t=o = kbT/m. The analytical solution for 
this case was found already in [13, [I4I and is given by 
(x^it)) = 2{kbT/m)t^ E2-a,3{-^t^~")- We find that for 
a < Uc {x'^{t)) is non-monotonic, while for a > ac it is 
monotonic. The non-monotonic behavior of {x'^(t)) for 
small a is a manifestation of a cage effect. Our finding 
that the same critical exponent ac describes both the 
phase diagram of the fractional Harmonic oscillator and 
also the fluctuations of a particle free of an external force 
field, emphasizes the generality of our results. Although 
not explored here in detail, the same ac is likely to de- 
scribe motion in any binding potential field, since it is 
the medium inducing the oscillations. We note that our 
results can be derived also from the fractional Kramers 
equation 1^, 22| and hence are not limited to the FLE. 

The response of a system to an oscillating time de- 
pendent field naturally leads to the phenomena of reso- 
nances, when the frequency of the external field matches 
a natural frequency of the system. The response of sub- 
diffusing systems to such time dependent fields was the 
subject of intensive research [2^, 24, [2^. In particular 
fractional approach to sub-diffusion naturally leads to 
anomalous response functions commonly found in many 
systems e.g. the Cole-Cole relaxation [2^. Here we 
investigate the response of a particle in a harmonic force 
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field, as found in the experiments 0, Hi, to a time de- 
pendent force Fq cos(r2t) using Eq. In the long time 
limit, using the stationarity of the process we find 



(x) — ^i?(r2)cos [nt + 9{n)] 

m 



t — > oo, 

where the response function is 



(16) 



(17) 



which was derived already in |27| using a different model. 
As well known for the normal case a = 1 a resonance is 
found when the friction is weak, namely a peak in R{U) 
is observed when u) > 7/\/2. On the other hand if the 
damping is strong uj < "f/V^ the peak is not found. The 
simple phase diagram describing the response of the frac- 
tional system is shown in Fig. [3l where a resonance phase 
means a phase with a peak in the response function. For 
a smaller than a critical exponent an — 0.441... one al- 
ways detects a resonance, even if a; — > 0. This is clearly a 
manifestation of the cage effect and it clearly shows that 
the critical exponent an marks a dynamical transition of 
the response of the system to an external force field. 

Many works consider the over-damped approximation 
to anomalous diffusion, which means that Newton's ac- 
celeration term is neglected, m = in Eq. ([1]). This ap- 
proximation must be used with care. As we showed when 
a < ac the exact solution of the FLE always exhibits a 
non-monotonic decay of Cx{t), while the over-damped 
approximation gives a monotonic decay. And according 
to our results when a < a^ resonances always appear in 
the response of the system to an external oscillating field, 
which are not found with the over-damped approxima- 
tion. In this sense the over-damped approximation fails 
when a is smaller than the critical exponents of the sys- 
tem. This is in complete contrast to the usual Langevin 
equation with white noise, where the over-damped ap- 
proximation captures the main features of the dynam- 
ics. This striking difference between the FLE and the 
usual Langevin equation stems from the fact that for the 
FLE the dissipative memory induces oscillations, as men- 
tioned. 

In conclusion, critical exponents mark sharp transi- 
tions in the behaviors of systems with fractional dynam- 
ics. The critical exponents describe a wide range of phys- 
ical behaviors: the correlation function Cx{t) of a particle 
bounded by a Harmonic field, the mean square displace- 
ment of the free particle, and the response of the system 
to an external oscillating field. Thus these critical ex- 
ponents are clearly very important and general in the 
description of the anomalous kinetics. The phase dia- 
grams we obtained are related to a cage effect, where for 
small enough a the medium induces oscillations in the 
dynamics of the particle. Thus the fractional dynamics 
is profoundly different from the dynamics described by 
the Langevin equation with white noise. 



RESONANCE PHASE 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0. 



FIG. 3; Phase diagram of the response of the system to a 
oscillating time dependent force field. Two simple behaviors 
are found either a resonance exists or not. For a < aji = 
0.441... a resonance exists for any binding harmonic field and 
any friction. 
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